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Abstract 

Based on a simple observation that a classical second order differential equation may be 
decomposed into a set of two first order equations, we introduce a Hamiltonian framework 
to quantize the damped systems. In particular, we analyze the system of a linear damped 
harmonic oscillator and demonstrate that the time evolution of the Schrodinger equation 
is unambiguously determined. 
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The search for field theoretic formulations of fundamental anyons 1 ' 2 has refocused at- 
tention on the long outstanding problem of the quantization of (non-Hamiltonian) systems 1- 
in which only the equation of motion is explicitly known. The common feature of this class 
of systems is its history dependence. The most notable such system is certainly that of 
the damped harmonic oscillator. Apart from requiring the knowledge of quantum aspects 
of this class of systems for practical reasons (for example, the laser system), this problem 
has its own right to be investigated both for the mathematical interests and exploring the 
new (perhaps deeper) fundamental features of the quantum theory. 

There have been quite a number of attempts to solve this problem, which may be 
classified into two classes. One is to start with a modified quantization scheme 4 ' 5 or a bigger 
system with introduction of additional degrees of freedom 6 such that the conventional 
Schrodinger or Heisenberg dynamical description is valid. The other 7 is to formulate the 
quantum theory in a modified Schrodinger dynamics (a nonlinear one). Clearly, attempts 
in the second class do not respect the superposition principle as the nonlinearity enters as a 
consequence of wavefunction dependent potentials. Similarly, the first class of the attempts 
also has some disadvantages. The introduction of the non-hermitian or time-dependent 
Hamiltonian in this class often leads to the conclusion of non-existence of the Schrodinger 
wave description. 4 

As is usually done in quantum mechanics, the time evolution of the Schrodinger dy- 
namics may be carried out only if one insures the completeness of the Hamiltonian eigen- 
states. The hermitian Hamiltonian in our standard theory guarantees unique unitary time 
evolution. For dissipative systems, the time evolution is certainly no longer unitary, yet one 
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may hope that the Schrodinger quantum mechanics may be still used for unambiguously 
determining the time evolution. The purpose of this paper is to present an alternative 
way of modifying the quantization framework in which the Lagrangians are introduced 
naturally and to show that the Schrodinger dynamics is indeed achieved. Thus the above 
mentioned puzzle will be resolved. # 

We shall begin with a general damped system on the line in an external potential, 
whose classical motion is govern by the non-linear second order differential equation: 

x + k(x)x + g(x) = 0, (1) 

where the dots denote time derivatives and k(x) and g(x) are real functions. We observe 
that a reduction of Eq. (1) to 

x + f(x) = (2) 

may be defined if / satisfies 

f(f'-k)+g = 0, (3) 

where the prime represents an x-derivative. Note the complex conjugate /* of / will be 
also a solution of Eq. (3) if / is a solution. One may verify this reduction by directly 
differentiating Eq. (2) with respect to time and using of the condition (3). Roughly 
speaking, this reduction reflects the fact that a second order differential equations may be 
decomposed into a set of two first order equations. For the obvious reasons, this reduction 
# During the preparation of this paper, we receive a paper by Feshbach and Tikochinsky, 5 
in which this puzzle was also resolved in a different way. We thank Roman Jackiw for 
bringing this reference to our attention. 
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has at least two advantages. Classically, a first order equation (linear or not) is of course 
simpler than the second order one,* while quantum mechanically, the first order equations 
may be regarded as the Hamiltonian equations by properly defining the canonical variables. 
For later convenience, we denote Xi(i = 1, 2) as the two (independent) solutions of Eq. (1) 
or the solutions of the follow set of the first order equations: 

Xi + fi(xi) = 0, (4) 

where fi are two (independent) solutions of Eq. (3). One may easily show that x = x\ +X2 
will be a solution of Eq. (1) only for the linear damped case, namely, /c=constant and g is 
a linear function of x. 
The Lagrangians L 

2 

L = J^iP* _ x *Pi + 2 fiPi]> ( 5 ) 
i=l 

which may be naturally introduced and give rise to Eq. (3) and 

Pi ~ ^{f'i{ x i)Pi + Pifl( x i)) = 0, (6) 

where no summation is assumed for the repeat index and the symmetric form is used 
for the usual reason. Thus the Hamiltonian is obtained in the usual manner (we use the 
symmetric form to deal with the ambiguous quantum ordering): 

H = -J2(fiPi+Pifi)- (7) 

i 

The symplectic two-form for the system is 

u> = 2 dpi A dxi, (8) 

i 

* These ideas may be used to solve the sine-Gordon and other nonlinear equations. 8 
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which defines the Poisson brackets. Therefore, the commutation relations are 

[xi,Xj] = \pi,Pj] = 

Using the rules (9) and the Hamiltonian (7), one may easily verify that the classical 
equations (4) and (6) are quantum mechanically realized (up to the quantum corrections 
0(h)) as the Heisenberg equations: 

x j = 7T l x j ' H] , 

\ (10) 

Pj = ^ H \- 

Note that Xi and pi may not be considered as hermitian operators and in the case of fi = f% 
the consistency condition requires the identifications x\ = x* 2 and p\ = p\ classically (in 
quantum case, they are x\ = x\ and p\ — p\). We will keep this in mind for later discussion. 

Now let us analyze the case of a linear damped harmonic oscillator, namely k = 
2A > and g(x) = cu 2 a;(mass=l). Thus we have / = — r\x from Eq. (3), where rji = 
— A± VA 2 — cu 2 [?7i(2) with the upper (lower) sign] with the convention r]i > 1]2 (or Im i]i > 
Im 772), if Vi are rea l ( or complex). Namely, the square root in this paper will be always 



positive and \/A 2 — u 2 = i^/u 2 — A 2 for A < u. The equations of motion (4) and (6) then 
become 

Xi T]iXi — 0, 

(11) 

Pi + ViPi = 0. 

The p-system is clearly related to the x-system by a time-reversal transformation U t : t — > 
—t, or equivalently U\ : A — > —A. Note that since the complex index % or the square root 
is introduced in our reduction, the time reversal transformation has been identified as the 
complex conjugation or simply shifting to in the above. Therefore the combined 

system (x = x\ + X2,p = Pi + P2) may be viewed to be time-reversal invariant. 
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The Hamiltonian H is from Eq. (7) 

H = ih(i] 1 N 1 + r] 2 N 2 ), 

which is indeed invariant under U a (a = t, A) and Ni(i = 1,2) are defined as follows 
Nt = -(atbi + bttti) = - X -{A % B % + BiAi) 

di = iAi, hi = iBi, Xi = a i exp(?7 i t), p { = ^ihbi exp(-r}it) . 

One may easily verify the following commutation relations using the rules 

[a,i, a,j] = [bi, bj] = 0, [ai, bj] = 5ij 

[Ai, Aj] = [Bi, Bj] = 0, [B i ,A j ] = 8 ij 

[Ni, dj] = -diSij, [Ni, bj] = bi5ij 

[N^ Aj] = -AiStj, Bj] = Brfij 

[N U N 2 ] = 0. 
From the relations (14), we immediately have 

>= ±K (±) + J)K (±) > 



«,K (+) >= V^P|nS +) -l>, 



h\n\ +) >= Jn\ +) + l|n (+) + 1 >, 



Bi\n\ ) >=Jn ( i Vi } -!> 5 



Ai\n\ ] >= ynj } + l\n\ } + 1 >, 
<nf\n'i >=5 n ^ i 6 cea/ (« = +,-) 

(n 4 (±) = 0,l,..-). 

Therefore the Hamiltonian H is exactly of the harmonic oscillator type and has the 
wave functions (of course, with the completeness) . 
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Note that since there is no prior reason to identify (bi) as the annihilation (creation) 
operators, we have used = iai and Bi = ibi (these transform the positions into the 
momenta and vice versa) to obtain the second set of solutions. Namely Bi (Ai) are our new 
annihilation (creation) operators. Alternatively, we may simply identify A\ = —ai as new 
creation operators. All these are similar to the usual case in which the canonical momentum 
operator consists of the operations of both annihilation and creation. Futhermore, it is 
worth pointing out that we have a similar situation to Ref. [5]. An operator, typically 
A = xp+px, [x,p] = ih (in our case, XiPi+piXi), seems to be hermitian but posses the pure 
imaginary eigenvalues. We resolve this apparent contradiction in a similar way to that in 
Ref. [5], namely, giving up the hermitian condition of A. Actually, an alternative method 
may also be employed, in which one may directly solve (—2ixd x — i)hip = Eip and assume 
ip to be single-value. Immediately, one gets E^ = ±i(2n ± + l)/i(n ± > 0, integer). The 
completeness of the states will be followed as the coherent type weight functions appear 
in the measure of the inner products. 

The above discussion also agrees with the representation theory as below. If one 
introduces 



one may easily verify that Jp(f3 = 1, 2, 3) are generators of the three dimensional Lorentz 
group SO(2, 1) or equivalently the two dimensional symplectic group Sp(2, R): 




(16) 



(17) 
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According to Bargmann, 9 all irreducible unitary representations of S0(2, 1) can be clas- 
sified as follows: (/j, < 0), (/i < 0), principal series, supplementary series, and 
the identity representation, where \i is the Casimir parameter, i.e. the Casimir operator 
jf — J\ — Jf = + 1). In terms of the A's and S's, there will an additional minus 
sign in the J's. Because of the sign change of the Jg, — J3 = \{AiBi + BiAi) will have 
the bounded from below spectrum as the square integrable wave functions are chosen. 
Thus, AT(±) = ±(n + \) as (i = =F^ (p, is determined by the commutation relations of 
a's and b's or A's and B's.), respectively. Note that for other values of fj,, one obtains 
the representations for one dimensional anyons, 10 or for parabosons, 11 depending on one's 
viewpoint. 

Clearly, the spectrum of H contains all information of the damped system Xi and its 
time reversed partner (growing one) Pi, as one expects. However, the spectrum of H is lost 
the usual meaning of energy levels and H may only be viewed as the dynamical generator. 
To be more specific, we focus on the case of the underdamping, A < u. In this case, the 

spectrum of H can be explicitly displayed as 

Ei+l = hn(n { 2 +) - n[ +) ) - th\(n[ +) + n { 2 +) + 1) 

(18) 

EtX = hn(n[ ) -n { 2 ) )+ ih\(n[ } + n { 2 } + 1), 
where O = \j u) 2 — A 2 and E^l 2 describe decaying states while En~^l 2 describe growing 
states. One may easily demonstrate that the states with E^ is indeed the time reverse 
of the states with E^~\ Our results (18) agree with those appeared in Ref. [5]. Since 
the well-known harmonic oscillator wave functions are complete, the time-evolution of the 
Schrodinger equation can be solved without trouble and of course is no longer unitary due 
to the damping force. For the harmonic oscillator, we see that the states of H becomes 
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those of the usual undamped harmonic oscillator in the limit A — > if only those states, 
satisfying the conditions ai\ip >= 62!^ >= are considered. This is similar to neglect 
the negative energy states which are physically unstable as we interpret H as the energy 
operator. 

Finally, we explicitly present here the Schrodinger time evolution under H for a given 
initial state |^(0) >= V (±) c ( ±) (± ) \n[ ±) n { 2 ±) >, 



n (±) 



E_i E (±) t . (±) (±) ( 19 ) 
c n (±) n (±)e R |ni ^ ; >, 



where E^n 2 are spectrum of if. Alternatively, one may easily write down the Schrodinger 
time evolution of an initial state \phi(0) >= V (+) CnX-2 W »4 > which decays with 
time 

m >= cW 2 e-^^ +) -"i +) )- A (^ +) +^ +)+1 )|nl +) n( +) > . (20) 

The time dependent operator expectation values can be obtained from either the 
equations of motion or from the standard Heisenberg representations. For example, 

Xi =< l(>(t)\Xi\l(>(t) > 

(21) 

=< v(o)Mt)|V(o) >=< ^(0)k|V(0) > e**, 
where by definition Xi = exp(iHt /h)cii exp(—iHt /Ti) (from the Heisenberg equation (10)). 
Note the definition of the "bra" < ip(t)\ is changed to the time reversed conjugation of 
\ip(t) >, rather than the usual hermitian conjugation of \ip(t) >. This is natural since the 
Hilbert space is defined by our time reversal invariant operator H. The same definition 
of the inner products appeared also in Ref. [5]. Likewise, one may easily write down the 
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correlations 

< iP(0)\ Xl (t)x 2 (t')\iP(0) >=< V(0)|aia 2 |^(0) > e" lt+ "< (22) 

In conclusion, we have analyzed the linear damped harmonic oscillator using the 
Hamiltonian framework, which is naturally introduced for describing this system and its 
time reversed partner. We have shown that the conventional Schrodinger dynamics is 
explicitly realized. Furthermore, the Schrodinger wave description has the usual interpre- 
tation and our method may be easily generalized to quantize three dimensional damped 
systems. Second quantization may also be carried out so that a field theory for this model 
is realized. 

Since our analysis starts directly by quantizing the classical solutions or phase space, 
as we have demonstrated, our method may yield a straightforward framework to study 
general non-Hamiltonian systems. We outline here how this may work. One will need first 
to identify the asymmetries which are responsible for the non-existence of the Hamiltonian 
of a system. These asymmetries are usually related to the time-reversal transformations 
because non-Hamiltonian systems have the common feature of history dependence. Then 
to find the classical solutions of both the systems will be the next step. Finally, one 
may find the Hamiltonian in which the original asymmetries will be restored to describe 
the combined system and the canonical quantization will be carried out directly. One 
such immediate example is the field theory of anyons. In 2+1 dimensions, the positive 
and negative spin corresponds to the different models 1 ' 2 and they are related by the time 
reversed symmetry (similar to the magnetic like interactions), 2 one may hope that our 
method may be used for quantizing anyonic field theories in a similar manner. 
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